Higher order contributions to Rashba and Dresselhaus effects 
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We have developed a method to systematicaUy compute the form of Rashba- and Dresselhaus-hke 
contributions to the spin Hamiltonian of heterostructures to an arbitrary order in the wavevector 
k. This is achieved by using the double group representations to construct general symmetry- 
allowed Hamiltonians with full spin-orbit effects within the tight-binding formalism. We have com- 
puted full-zone spin Hamiltonians for [001]-, [110]- and [lll]-grown zinc blende heterostructures 
{D2d,CAv,C2v,Cz-u point group symmetries), which are commonly used in spintronics. After an 
expansion of the Hamiltonian up to third order in k, we are able to obtain additional terms not 
found previously. The present method also provides the matrix elements for bulk zinc blendes (Td) 
in the anion/cation and effective bond orbital model (EBOM) basis sets with full spin-orbit effects. 

PACS numbers: 73.21.-b, 71.15.-m, 71.15.Ap 



I. INTRODUCTION 

There have been strong recent interests in spin ef- 
fects in semiconductors for applications in spintronics^. 
Particularly, spin splittings due to bulk inversion asym- 
metrji^ (BIA) — arising from the different chemical char- 
acter of the anion and the cation in the zinc blende 
structure — and structural inversion asymmetry^ (SIA) — 
appearing in layered structures where "top" is different 
from "down" : asymmetric composition, doping, etc. — , 
are studied because a good understanding of the intra- 
band spin splittings is essential in any attempt to under- 
stand the operation of spintronic devices at a microscopic 
level, since they determine the spin dynamics during re- 
laxation and transport processes^iSiSi^iifiiiiiiSiii. The 
pioneering work by Dresselhaus^, Bychkov and Rashba^r, 
and D'yakonov and Kachorovskii^ taught us the func- 
tional form of the leading order contributions — up to 
first (third) order in the wavevector k for SIA (BIA) — to 
the spin Hamiltonians in bulk and heterostructure zinc 
blendes due to the various sources of inversion asymme- 
try. 

This constructive procedure is of great use as it de- 
scribes the main physics but, as shown below, 0{k^) con- 
tributions for heterostructures due to BIA^ should con- 
tain additional terms, while the 0{k^) contributions aris- 
ing from SIA have not been studied. Thus, it is also im- 
portant to have a systematic modeling tool which guaran- 
tees that all spin-related qualitative features in the band 
structures will be present. The empirical tight-binding 
method as formulated by Slater and Koster— provides 
such a systematic way of generating all the symmetry- 
allowed terms that can appear in a Hamiltonian. It 
has been used extensively in the computation of bulk, 
heterostructure and surface properties, and today finds 
widespread use in the study of nanostructurcs involv- 
ing thousands and even millions of atoms^^, or transport 
properties'^. 



However, the original formulation by Slater and Koster 
was obtained through the use of single group symmetry 
operations, which prevented it from describing spin ef- 
fects. Later, Chadiii extended the original method with 
an on-site spin-orbit energy which effectively describes 
the zone center split-off splitting and breaks the dou- 
ble degeneracy of the bands for systems without an in- 
version center. This scheme contains the essential fea- 
tures of spin-orbit interaction, but does not reproduce 
some qualitative aspects such as the linear spin splitting 
in valence bands in zinc blendes^. Boykini^ has intro- 
duced supplementary matrix elements between nearest- 
neighbor atomic orbitals as a remedy for zinc blende 
structures, but that procedure does not guarantee a pri- 
ori that all spin-orbit effects will be considered, and does 
not treat systematically other types of structures. The 
development presented in this paper allows us to address 
these issues precisely. 

In this Article, we compute the forms of SIA and 
BIA contributions to the spin Hamiltonian of com- 
mon heterostructures, beyond the conventional Rashba™ 
and Dresselhaus^ terms, to an arbitrary order in the 
wavevector k. This allows us to treat SIA on the same 
footing as BIA, providing higher order corrections to 
the Rashba Hamiltonian for [001] {D2d,Ciy,C2v)i [HO] 
{C2v,Cs) and [111] (Ca^) diamond and zinc blende quan- 
tum wells. We achieve this by following the constructive 
process of Slater and Koster^"* (SK), but using double 
group irreducible representations (irreps) as opposed to 
single group irreps. We also apply the method to the gen- 
eration of bulk zinc blende [Td] Hamiltonians in the an- 
ion/cation and effective bond orbital model (EBOM)i2i2£ 
basis sets containing full spin-orbit effects. The proce- 
dure presented could, of course, also be used for the de- 
scription of the bands of materials involving heavy ele- 
ments (i.e. large relativistic effects), such as lead com- 
pounds, rare earths, etc. 
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II. METHODS 

In order to achieve our goal of constructing spin Hamil- 
tonians to arbitrary order in the wavevector k, we first 
construct the tight binding Hamiltonian for the corre- 
sponding point group symmetry using the EBOM basis, 
and then series-expand in k to the desired order. Also, 
since we are interested in the correct description of rel- 
ativistic and, in particular, spin effects in the electronic 
bands, we are led naturally to the use of double group 
irrepa^ to describe the symmetry operations of the crys- 
tal. Thus, if we follow the SK procedureM using double 
as opposed to single groups, we will be assured to ob- 
tain the most general Hamiltonian compatible with the 
crystal symmetries which, by construction, will include 
all spin-orbit effects. 

The basis states of double group irreps are much less 
amenable to brute force operation than their single group 
counterparts. Therefore, it becomes necessary to write 
the tight binding equations and the symmetry constraints 
of the parameters in a form that allows computer ma- 
nipulation. This is done in Appendices fXl and IBl which, 
while rederiving the empirical tight-binding basic results, 
also set forth the notation that will be employed through- 
out the rest of the paper, and help establish a better 
connection with previous authors' conventions. For our 
actual double group calculations we use the matrices 
D{G)1g as provided by Shirai^^, which correspond to the 
ones tabulated by Onodera and Okazaki2i. 



III. RESULTS 

In order to illustrate the procedure with a familiar 
example, and for comparison purposes, we construct in 
Sec. nil Al a 4-band model with no spin included for a 
zinc blende (point group Td) using an effective bond or- 
bital basis set. Sec. IIII Bl presents the main results of our 
work, with spin Hamiltonians for various sorts of quan- 
tum wells displaying the terms additional to the Rashba 
and Dresselhaus Hamiltonians. Finally, we show explic- 
itly in Sec. IIII Cl the equivalence between the Lowdin or- 
bitals (LOs) centered at the primitive cell sites [effective 
bond orbitals (EBOs)] and the zone center solutions of 
the k • p theory. 

The definition of the symbols appearing in the follow- 
ing tables can be found in Equations ljA4|l and (jA5|) in 
Appendix IXI The phase of the parameters Eij{r) — the 
Hamiltonian matrix elements between a state i at the 
origin and a state j at r — has been factored out in the 
calculations of the tables. Thus, all Eij (r) are real. Pa- 
rameters appearing in different tables are of course un- 
related, while the procedure in Appendix ^ guarantees 
a high degree of independence within the parameters in 
any one table. 



//f/ri(k) £,,(000) + 

4Ess ( 1 5 0) (cos ^ cos 7] + COS r/ cos ( + cos ^ cos ^) 
H^i (k) 4iS,, ( i i 0) (cos -h cos C) sin C - 

4Eszi^^0) sin sin C 
H^f^- (k) E.^^ (000) + AE^^ ( i iO) (cos n + cos C) cos ^ + 

4£^^(ii0) cos 7? cos C 
Hlf^ik) -4i5,,(ii0)sinCsin77 + 

4iSa;z(ii0) (cos^ — cos 77) sm( 



TABLE I: Matrix elements for the zinc blende structure (sin- 
gle group). The definitions ^ = kxa/2, r] = kya/2, = kza/2, 
are made. The remaining symbols are defined in Equations 
llA4ll and llASI l. For example, Hg^^^{]s.) represents the ele- 
ment of the Hamiltonian connecting a state belonging to the 
Fi irrep with s symmetry to a state belonging to the F5 irrep 
with X symmetry. The parameters Eij(r) appearing in a ta- 
ble have no relationship with similarly named parameters of 
a difi'erent table. 



A. Single Group 

We will study the top valence and lowest conduction 
bands of a zinc blende (point group T^), which have F5 
and Fi symmetry at the zone center (we will use the 
KDWS notation^^ throughout this paper for the irrep 
labels). We consider only coupling to the twelve nearest- 
neighbor sites. Since we are working with effective bond 
orbitals located at fee lattice sites, the nearest neigh- 
bors are at (55O) and its equivalent positions. Table U 
shows the matrix elements. The use of the arguments in 
Appendix IbI shows for example that, with our choice of 
phases, the whole block E^^^'^ [see Eq. (jA5|l for mean- 
ing] is purely real — which is otherwise trivially obtain- 
able since the Hamiltonian without SO is real and the 
p orbitals have the same (imaginary) phase — and that 
EzxiW^) ^ —ExziW^)- These results coincide with 
those of Hass et ali^, which correct the misprints in 
Table V of Ref. or those obtained by the different 
method of adding a d component to p states in an fee 
lattices^. 



B. Double Group 

In what follows, we present models of quantum wells 
(QWs) of different symmetries, including structural in- 
version asymmetry (SIA)^ effects only, bulk inversion 
asymmetry (BIA)^ only, or both; having in mind the 
study of in-plane spin transport in two-dimensional elec- 
tron gases (2DEGs). This has the advantage of yielding 
less cumbersome expressions than those of the superlat- 
tice of the same symmetry, which can be obtained if nec- 
essary by setting the matrix elements between supercell 
instances of different z, Ef^ [l, m,n^ 0), to a finite value. 
Therefore there will be no terms in the Hamiltonian 
matrix elements. 

Of course, for a numerical zinc blende nanostructure 
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computation the starting point would be a Hamiltonian 
derived from the bulk matrix elements given at the end of 
this section. In this context, the following results should 
rather be taken as an analytical tool for the study of the 
bands and the spins, and they can be used to test whether 
calculations from other methods satisfy the symmetry 
requirements. 



H^f" £'tt(000) +2£tT(100) (cosf + COS77) + 

4£'||(110) cos ^ cos 77 
H^eTe -2En{W0) (sin^ + isinr?) + 

2-/2-E| i ( 1 10) (sin ^ cos 77 + i cos ^ sin 7]) 



TABLE III: As Table El but with BIA-only effects (L>2d). 



1. [001] quantum wells 

a. SIA only In this configuration only the Rashba 
splitting^ should appear. These quantum wells (QWs) 
possess Civ symmetry, and an example could be a 
Sii_a;Gea;/Si QW with an asymmetric Ge concentration 
profile. For C^y QWs the CB will transform according 
to Tq. We construct a model for the CB with on-site and 
second nearest neighbor coupling, with the assumption 
that the supercell instances do not couple to each other, 
and the results are shown in Table UTI 

Expanding again about the F point and keeping terms 
of up to third order, we obtain 



spin 



2a 



i?Ti(100) + V2i?t4(110) {ky (7x kxCy) H~ 
£;tx(IOO) + V2E^i{110)\ aV3 (-fc^ffx + k^ay) + 

V2a^ETi{U0)kj,ky{-k^ax + ky<Ty). (I) 



The first term of the right-hand side is, of course, the 
well-known Rashba Hamiltonian. Here we also derived 
the two 0{k^) terms, which are the next order in impor- 
tance. Up to present there have been calculations that 
have been carried out to third order in BIA effects but 
only to linear order in SIA terms'". The results from 
these calculations can be improved and be made con- 
sistent within the chosen order of approximation if the 
higher order SIA terms are included. 

b. BIA only Here the corresponding point group for 
a zinc blende [001] QW is D2d- The different orientation 
of the bonds at the interfaces [native inversion asymme- 
try (NIA)22i] may lower the symmetry to C2V, but this 
case is equivalent to SIA+BIA, and is treated in the next 
section. Applying Eq. IjAlOp to the CB (Fg) we obtain 



I (000) +2E„{ 100) (cos C + cos 77) + 
4iJ||(110) cos ^ cos 77 
2£;ti(100) (j sin 5 sin 7?) + 
2^/2S||(110) {i sin^ cos rj + cos 4 sin 77) 



TABLE II: Matrix elements for an [001] structure with SIA 
only described by the dv point group. The definitions ^ = 



_Bt|(000) + 2£n(100) cos^ + 2^11(010) cosr; + 
4i5||(110) cos 5 cost; 
H^f"" 2£;||(100)isinC - 2£;||(010) sinr; - 
43 [£;ti(110)]cos^sin77-(- 
4K [E^ ^ {noy\ i sin ^ cos ri 



TABLE IV: As Table lul but with SIA and BIA effects {C2v). 
The phase of i5||(110) has not been factored out as it is not 
determined by symmetry. 



the results in Table UTTl and expanding close to F: 



spin 



{-kxUx + kyGy) 

\/2a^ Ei^ l{ll{))kxky {-kya^ + k^cFy) + 
£;Ti(100) - V2En{llG)\ aV3 {kla^ - k^ay) . (2) 



We recover now the Dresselhaus Hamiltonian for [001] 
zinc blende quantum wellsSiSSiSS in the first two terms. 
The middle term can be obtained by the procedure of 
taking the spin Hamiltonian for bulk zinc blendes [Ref.H 
and Eq. (|13|l below] and substituting fc^ for its expecta- 
tion value in the quantum wcll^. However, the last term 
cannot be obtained in this fashion^SS which shows the 
power of our proposed approach. The symbols i?|-f(100), 
i?|x(110), etc., have no relationship with their analogous 
in the other cases, and the reference to the corresponding 
point group has been dropped to lighten the notation. 

c. SIA+BIA The corresponding point group for 
these structures is C2V , and the conduction band is asso- 
ciated with the F5 irrep. Table Hvl shows the results, and 
the spin Hamiltonian is 



i7. 



kxU and 77 ; 



kyO are made. 



« ^2a{[En{0l0) + 25 [EniUO]]] kyax+ 

[i;n(ioo) + 25R [En{im]kxCTy} + 

aV3 { [En (010) + 23 [E^ (HO)]] klax + 
[i?Tx(100) + 25ft[£;u(110)]]fc>,+ 
63 [STi(llO)] klkya^ + 63? [^Ti(llO)] k^k^ay} . (3) 



This Hamiltonian can be transformed into a more fa- 
miliar form by making the substitutions En(OlO) — )■ 
Ebia + EsiA, i^Ti(lOO) ^ £;bia - ^^sia and £^11 (HO) ^ 
(£^BiA,iio + -EsiA.iio) /2 + i (i?BiA,iio " -^SiA.iio) /2 and 
reverting to the original axes definitions for the zinc 



FIG. 1: Labeling of the irreducible representations (irreps) corresponding to the lowest conduction band and the two highest 
valence bands, at the zone center, for structures with different point group symmetries. Spin splittings at A: 7^ are omitted 
for clarity. 



blende structure: 

Hspil" ~ 2a [£;siA - 2£;siA,iio] {-kya^ + k^cTy) + 
2a [Ebia + 2i?BiA,iio] {-kxCTx + kyCTy) + 
aVe {[£;siA - 4£;siAaio] {kla^ - klcTy) + 
[Eb,\a + 4£'bia,iio] {kl 
SEsiAkxky {kxCTx ~ ky(Ty)+3EBiAkxky [kydx ~ kxffy)} , 

(4) 

which contains the same functional dependence as the 
sum of the separated SIA- and BIA-only cases, Eqs. 
and 

2. [110] quantum wells 

a. SIA only The symmetry group is C2V, so the re- 
sults are essentially the same as for [001] heterostructures 
with BIA+SIA. It is worth noting that, because of the re- 
duced symmetry with respect to SIA- [001], there appear 
additional Dresselhaus-like terms although for SIA-only 
one would expect only Rashba terms. 

We have verified these predictions with the help of k • p 
code where BIA effects can be turned on and off"^^ for 
a [110] 16/16 AlSb/GaSb/InAs/AlSb QW (see Fig.[2J. 
We see for this case that the splitting presents a small 
degree of anisotropy due to the interplay of Rashba- and 
Dresselhaus-like terms in the spin Hamiltonian, although 
the simulation only accounts for SIA effects. The axes 
of constructive and destructive Rashba-Dresselhaus in- 
terference are rotated by it /A with respect to the [001] 
BIA-I-SIA caseSS because here the planes of symmetry 
coincide with the usual choice of x and y axes, while in 
the [001] the reflection planes bisect the x and y axes. 

b. BIA only The symmetry group is again C'2v , but 
now the twofold axis is in-plane, whereas in the previous 
case it was along the growth direction. This changes the 
symmetry considerations considerably. The results for a 
second-nearest-neighbor model are displayed in Table M 

An interesting feature that is recovered is that, for — 
(or for a quantum well) , the spin will align or antialign 
along the growth direction. This result still holds even 
when all the neighbors are included: a general k point in 




270 

FIG. 2: Polar plot of the spin splitting in the CB for a [110] 
16/16 AlSb/GaSb/InAs/AlSb QW with SIA only calculated 
with k ■ p. The anisotropy of the splitting is due to the inter- 
play of Rashba- and Dresselhaus-like terms in the spin Hamil- 
tonian, although the simulation only accounts for SIA effects. 



H^f" -Ett(OOO) + [£n(100)] cos^ - 

23 [^1 1 ( 100)] sin f + 2£n (010) cos 77 + 

4K [E^^lnO)] cos ^ cos 77-49 [£||(110)] sin^cosr? 



TABLE V: Matrix elements for an [110] structure with BIA 
only described by the C2v point group. The definitions ^ = 
kxUx and r/ = kyUy are made. 



the well Brillouin zone has Cs symmetry (reflection with 
respect to the well plane), which confines the spin to lie 
perpendicular to the plane. This makes symmetric [110] 
structures have unusually long spin lifetimesS*^ for spins 
along the growth direction because of the suppression of 
the D'yakonov-Perel' (DP) spin relaxation mechanismA 
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H 



TT 



H 



n 



£;tt(ooo) + 

23? [£;||(010)] (2cos V3C/2cosr?/2 + cos77) + 
23[_Bn(010)l (2cosx/3C/2sm77/2 - sin??) + 
2£n(V300) [cos ^5 + cos(V3C/2 - 877/2) + 
cos(V3C/2 + 3?7/2)] + 
-2£'||(010)(cos V5C/2 sin 77/2 + sin77 + 
^^/3 sin ^^/2 cos 77/2) + 
iEuiVSOO) [2smV3^+ (1 + 
1^3) sin - 37?/2) + 

(1 - iVS) sin(V3?/2 + 37;/2)] 



TABLE VI: Matrix elements for an [111] structure described 
by the C^v point group. The definitions ^ = k^a and 77 = kyU 
are made, and the primitive vectors are taken to be a(0, 1) 
and a(^/2, 1/2). 



Then, the spin Hamiltonian takes the simple form 



H, 



« -2a,3 [E^^ (100) + 2E^^ (110)] k,a,+ 
a,/3 {al^ [i^TT (100) + 2E^^ (110)] kl+ 

0^63 [E^^{nO)]kl}k,a, 



(5) 



c. SIA +BIA In this case the symmetry will be low- 
ered to Cs (a single reflection plane), and the above argu- 
ments will not be able to put any restrictions on the spins. 
This is further understood if we combine the BIA-only 
results, which keep the spins pointing along the growth 
direction, with the SIA-only results, which tend to keep 
the spins perpendicular to the growth axis. Thus, the 
spins will point in a general direction dependent on the 
particular structure under study. 



3. [Ill] quantum wells 

For BIA-only, these quantum wells transform accord- 
ing to the Czv point group'^^; consistent of the iden- 
tity, two threefold rotations about the growth axis and 
three reflection planes separated by 120° that contain 
the threefold axis. It is easy to see by inspection that 
the inclusion of SIA does not reduce the symmetry with 
respect to the BIA-only case. Thus, here it is only nec- 
essary to derive the TB Hamiltonian for Cst, structures, 
which is shown in Table IVTI 

The expansion about the F point yields a spin Hamil- 
tonian: 



Ht. 



3a 



kxCFy) 



\/3i;|x(%/300) -£;||(010) (kya, 

a^/sj 3i;|x(010) - 9V3i;ti(V300) k^ {kyOa, ~ hj^ay) + 

23[£;n(010)](-3fc2 + /c2)fc,a,}. (6) 

Thus, no matter for what situation in [111] QWs, we 
only have — at 0{k) — Rashba-like terms, as shown previ- 
ously with a different methodic. This is important be- 
cause the Rashba contribution can be tuned through the 



action of a gate bias, and it is in principle possible to 
make the first order spin splitting vanish, suppressing 
the DP spin relaxation mechanism for all components at 
the same timei^. 



4.. Zinc blendes with effective bond orbitals 

Here we will construct the model equivalent to the 
^rsFs |-,}qj,]j Sec. nil A1 but with spin. Thus, it will in- 
clude a heavy- and light-hole Fg and a split-off F7 set of 
states. We provide the symmetry-constrained matrices 
E^'^^ , as they are the primary quantities in nanostruc- 
tures: 



prgrs/i 1 
^2 2 



0) 



£33 


UJ*E3i 


^Eri 


uE^^ \ 


-UJE31 


En 


u;*E,i 


^£31 




-ujEj^I 


En 


tO*E3i 




-iE^i 


—L0E31 


E33 I 



E^^^^a^o) 



E 



and 



T7TS 



-U)E; 



-iE^ 



3T 



E 



3T 



''Cl'*_£/3-|- 


iE-^] \ 




uj*Ei^ 




-£iT 1 ' 


^ iE^ 


ioE^^ J 


Ei-\ 


— LO"^ E\^ - 


ujE\t^ 


E\^ — 


E]] 


uj*E^l\ 


-ujEii 


E]] ) ' 



(7) 



(8) 



iE. 



UJ*E: 



3T 



3T 



(9) 



(10) 



^[;2> 



and t (i) refers 



where u; = e'^/'*, I = —I, I refers to 

to \4\J^) {{(j/^^^^))^ s-iid the phase factors are exphcitly 
shown so that all parameters Eij are real. Thus, using 
symmetry operations, the 36 initial matrix elements are 
reduced to 12 independent parameters. The matrix el- 
ements E^'^^ (Grj) for the remaining nearest-neighbor 
sites are obtained combining Eqs. ((7|)- ((Tn|l and (|A8|) . 

Now we can construct the bulk tight-binding matrix 
elements, which will be the foundation for further anal- 
ysis of the meaning of the parameters. Table IVIll shows 
the computed matrix elements from Eqs. (|A10I) and (|B8|I . 
The remaining matrix elements can be obtained from the 
hermiticity of the Hamiltonian in the Bloch sum repre- 
sentation. We can obtain information about the physical 
effects that are included in this model by expanding the 
matrix elements^^ about k = and comparing to k • p 
results^!. 

We focus first on the heavy hole (HH) and light hole 
(LH) bands, given by the F7 irrep. ETB implementations 
based on atomic or bond orbitals with the addition of on- 
site spin-orbit ^^'^^ yield an incorrect intraband splitting 
proportional to k^ along [110] , while Dresselhaus showed^ 
from symmetry arguments that it should be linear with 
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H^f^{-k) ^33(000) + [3£ii(ii0) +^33(^0)1 (cose COS T] + COS rj cos ^ + cos C cos e) + 3i?33 (5 10) COS e COS ri + 

^[SEniUO) + 2%/3i53i(|i0) - i?33(||0)] (cosC - cos,)) sine 

//[f^**(k) -£33(000) + [SBsadIO) +-Bii(iiO)] (cos^cosT^ + cosrycosC + cosCcosC) +3-Bll(iiO)cosCcos?7- 

^[i5lI(liO) - 2y3i53i(ii0) - 3£;33(M0)] (cosC - cos ,7) sinC 

H^fHk) - [\/2i53i(M0) + \/W2(i?ii(M0) + S33(530))] (smC + jsmr,)cosC + 

2A/2i53i(iiO) (sin^ cost; + i cos ^ sin 77) + 4E-^i{^^0) (sin^ — isinr;) sin^ 

Hlf^ (k) [_y2£3i (i iO) + 7372(^11 (i |0) + i533 (HO))] (cos C + COS 77) sin C + x/3(i5n (i iO) - 

-B33(||0)) (cos^ — cost;) cosC — 'iiEgi{^^0) sin^sinr; 

Hlf^ik) [_3£ij(li0)/V2 + ^i?3i(M0)+ ^33(550)/^ (sine-7sin77)cosC + 

2V2-B33(iiO) (— sin 5 cost; + icosf sin?]) 
Hff^«(k) _[s^j(iio)/V2 + x/6S3i(iiO)-3£;33(iiO)/\/2] (sin ^ + i sin r,) cos C + 

2A/2i5ij(iiO) (sin e cost; + i cos ^ sin 77) 

^ (^^ ' ' ~ ^ i^x 5 7 ) 5 ^ {kx 5 J ^2) — ^ (kx ; ^1/ 1 kz^ 

H^f'Ck) Stt(OOO) +4£tt(ii0) (cos? COST? + cos 77 cost;" + cos ("cos ^) + 2\/2i?||(ii0) (cos 5 — COS 77) am( 

H^f^k) 2^_Bu(iiO) [(cos 77 - cosC) sin^ - i (cosC - cos^) sin 77] 

//Jf^^(k) 2V2£;3t(|50) (sin 5 cos 77 + 7 cos C sin 77) - {V2E3^{^^0) + a/G-Bi^ (i ifl)) (sin ^ + i sin 77) cos C + 

2£'3|(ii0) (sin^ — isin77) sin(" 

H^f^'' (k) 4_Bi| (i iO) cos ^ cos 77 - 2£;i| (^ ^0) (cos ^ + cos 77) cos C - (VGEa-ii^^O) + \/2£i| (i |0)) (cos ^ - cos 77) sin ( 

H^f' (k) 2y2£;i^ ( i iQ) (- sin e cos 77 + i cos ^ sin 7?) + (\/6S3t ( " ^-^iT ( M^)) (^in C - » sin 77) cos ( + 

2y3£;3|(ii0) (sin^ + isin77)sinC 
H^fHk) -4iB3|(ii0)sin^sin77-2V^Sil(ii0) (cos C - COS77) cos (+(^2-^31 (i i0)-^/6£i| (| §0)) (cos ^ + cos 77) sin C 

7^21* ^ (^a; 5 5 kz^ ~ -^31* ^ (^a: ; ky . ) ? ^ (^^ ' 1^2) ~ ^ (^^ ' ' ^2)5 

^11 ^ (^^ ' ' ~ ^IT^ ^ ' ' ) ' ^ (^^ ' 1 ^2 ) = ^ (^^ ' ' ^2)5 



TABLE VII: Matrix elements for the zinc blende structure (double group). The remaining matrix elements can be obtained 
from the hermiticity of the Hamiltonian in the Bloch sum representation. The definitions ^ = kxa/2, 77 = kya/2, (" = kza/2, 
are made. 



k. Our current method solves this problem. From k • p 
theory we kno"w that a parameter C describes the linear 



splittingii^. Looking, for example, at the H^^^^ matrix 
element and comparing to the corresponding k • p matrix 
element, ~^/3C{kx + iky)/2, we can identify 

C = V2a [VSEni^^O) - 2i;3i(ii0) + VSE^-^i^^O) 

(11) 

If the parameters in C were calculated from single group 
theory, we would obtain Eii{^^0) = -2V2E^,{^^0)/3, 
Esii^^O) = -v/273£;,,(iiO) and E^si^^O) = 0, yield- 
ing C — and, therefore, no linear splitting in the HH- 
LH bands. 

Figure 01 shows the calculated HH, LH and SO bands 
for bulk GaSb. For actual calculations to show the inclu- 
sion of the effect, we set the double group parameters to 
the single group values as determined from the Effective 
Bond Orbital Model (EBOM) methodic, with the addi- 

We have taken C = 0.7 



tion of £;33(iiO) = 
mcV-A as calculated by Cardona et al^ Plot (a) shows 
the full zone bands. The inclusion of £'33(^^0) makes no 
appreciable change to the bands at the scale of the plot 




b) 




LAP A X U,K 



-0.15 -0.10 -0.05 0.00 0.05 0.10 



k[10 2it/a] 



FIG. 3: a) Full zone HH, LH and SO bands for GaSb. b) 
Bands very close to the zone center along [100] and [110] di- 
rections calculated with the present method (solid lines), fea- 
turing the linear splitting, and calculated with the method of 
Ref. (dashed lines). 



{i.e. no linear splitting is visible at this scale). Plot (b) 
is a closer look at the HH and LH bands at the top of the 
valence band comparing the method presented here (solid 
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lines) and a previous EBOM implementation allowing 
for bulk inversion asymmetrji2& (dashed lines). Only the 
present method reproduces the linear splittings, in accor- 
dance with the predictions from group theory"^, showing 
that the linear splitting arises from off-site SO interac- 
tions, which were not included in R,ef.l2(il In this particu- 
lar case, the scale of the energy is tiny, and it makes little 
sense to make any claim when such small energy scales 
are involved. Nevertheless, this example shows that it is 
indeed possible to describe correctly all spin-orbit effects 
within a tight-binding framework and, moreover, the ef- 
fects of a finite C have been observed experimentally for 
GaSb in hole transport measurementsS&Si. We show be- 
low more examples (see, for example, Sec. lIIIB2l a) with 
potentially bigger magnitudes. 

We now turn our attention to the lowest conduction 
band in zinc blendes, which is associated to the Fg irrep. 
All the symmetry-equivalent points of (^^0) are spanned 
by the proper rotations of T^;, which under Fg correspond 
to the same matrices as T^. Therefore, the H^''^'^ matrix 
elements remain unchanged under the substitution Fy — s- 
Fg, leading to a full zone spin Hamiltonian 



with 7c = a^/{2V2)E^f''{^^0) and a, the Pauh spin 



matrices. 



5. Zinc blendes with quasi-atomic orbitals 

Within ETB, a more common approach for the study 
of zinc blendes is the use of orbitals assigned to atomic 
position3i^i2S4^2ii£ rather than Bravais lattice sites. This 
increases the spatial resolution and the size of the basis 
set. The anion is assumed to have da — 0, while the 
cation is at dc = a(l, 1, l)/4. The blocks for the first- 
neighbor overlaps for VB states are then given by 



E 



/T'r^riM 1 In 

^ V4 4 4^ 



fill) = 

U 4 4'' 











iV2E^3 --ujV3E^3 



-uj*V3E^3 iV2E^3 
uj*E^3 



(14) 



(15) 



H. 



2V2E^f''{^^0) [a^ (cosC - cos??) sinC + c.p.' 



(12) 

where c.p. stands for cyclic permutations and which, after 
expanding to lowest order about the zone center, allows 
us to recover the Hamiltonian^ 



H. 



7c [{kl - ky) k;,az + c.p. 



(13) 



^ U44/ 



and 



f Lu*E3^ -iV2E3^ ^^ 
-Lj*V3E3^ 
—cj^/SEj-^ 



(16) 



^ U 4 4^ 



£-33 

-oj[V2E3i+E3s/V3] 

-Z^31- 21^273^33 



-LU*[V2E3i^ 

E33 
-ujE. 



V^E; 



33J 



33 



-iE. 



31 



W*£'33 
-£'33 

Lu[V2E3j + ViE^s] 




(17) 



and E^i^i (^^^) *^^ii obtained from 

"'"£;,';;^(-R, - -I- d„) - ""'£;,^f'*(R, + - d„), 

(18) 

with a, a' = a,c (anion or cation). The spinless model is 
recovered setting i?33 = E-i^ = Exx, i?3i = E^] = E-^3 = 
Exy/V^ and E33 = 0. The non-trivial Hamiltonian ma- 
trix elements in the bulk orbital representation are given 
in Table EnH 

C. Effective bond orbitals and k • p 

We can obtain further insight into the nature of the 
Lowdin orbitals (LOs) used as a basis set by noting that, 



for orbitals associated with the primitive cells [effective 
bond orbitals (EBOs)] and at k = 0, we have 

E 1^ E ^(G)r,i?ri;'(r,)i?(G-i)^;, (19) 

and by use of the orthogonality theorem for the irre- 
ducible representations^ we arrive at 

il£f'(k = 0) = S,,,Su' Y: ^^^st'ir,), (20) 

where N{Oq) is the number of operations of the point 
group and is the dimensionality of the /x-th irrep. 
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//[f ^(k) i [s3i(iH) + VW2S33(H3)] (9.+i9y) 

<f«(k) iS3l(U^)5. 

<f«(k) i£33(lil)(g,_jgJ/V2 

H[3^^^(k) _i[i?3j(lii) + ^£33(111)] {g^~^gy) 

//[f^(k) J533(iii)5o 

H['f^(k) -ii533(iii)(5,+^)/V2 

<f «(k) i [i53i(iil) +2v'273i533(|H)] ff^ 
-'^ij* ^ i^^j ^vt ~ H^^^ ^ (kx, ~ky, ~kz) 

H^7^^(k) i5TT(iH)ffo 

<rj(k) 

papc papc 

H^f^^k) 

Jf[f^(k) _jvW2i53T(i3 3)(ff.-jff.) 

H[«^J(k) _iV2i53T(||i)ff. 

^il (^2; , fcj/, fcz) = — -ffjj* {kx, —ky, —kz) 

<f^(k) 

Jf[f^(k) j^£;^3(lll)(g,+igj 

H[|^J(k) j^£;^3(i|ijg. 



TABLE VIII: Matrix elements for the zinc blende struc- 
ture (double group) for quasi-atomic orbitals. The remain- 
ing matrix elements can be obtained from the hermiticity 
of the Hamiltonian in the Bloch sum representation. The 
definitions go = 4(cos ^ cos cos — isin^sinr;sin^), gx ~ 
4(— cos ^ sin 77 sin -)- i sin ^ cos r;cos C^) (and cyclic permuta- 
tions), are made, where ^ = kxa/4, rj = kya/4 and ( = k^a/A. 

Thus, if one is constructing a model where there are no 
repeated irreps, we see that the Bloch sums of EBOs will 
exactly diagonalize the Hamiltonian at k = 0. In other 
words, the Bloch sums of EBOs will be the zone center 
solutions, making the connection between the EBOs and 



the k • p basis evident. 



IV. SUMMARY 

In conclusion, we have constructed full-zone spin 
Hamiltonians for [001], [110] and [111] zinc blende quan- 
tum wells. We then performed small k expansions of 
those Hamiltonians about the zone center, yielding their 
k • p counterparts. The k • p Hamiltonians thus obtained 
present spin-dependent terms that had not been previ- 
ously described in the literature. In particular, we see 
that the Rashba Hamiltonian can be supplemented with 
third order terms, which will need to be included in cal- 
culations where other sources of spin splitting are con- 
sidered up to that order. We also generate additional, 
growth direction-dependent 0{k?) contributions to the 
Dresselhaus Hamiltonian. The method we have employed 
is not restricted to these particular cases, as it extends 
the tight binding formalism to include the treatment of 
spin by using double group representations. This guar- 
antees the systematic inclusion of all spin-related effects 
in the bands. Thus, our work can serve as the basis 
for numerical studies of large scale nanostructures where 
spin effects are important, as well as an analytic tool for 
predicting spin properties in reduced-symmetry systems. 
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where S~^^^ is the inverse of the square 
the overlap matrix S, with elements 



(Al) 



root of 



(<y9,^ ;Rj'|(pf ;Rj), and ;Rj) is the atomic orbital (or, 



in general, a state with finite overlap, e.g. a bond or- 
bital) analogous to |0f;Rj). This is well defined, as the 
Lowdin symmetrization procedure ensures that the state 
will have the same symmetry properties as the underlying 
atomic orbital from which it is constructed^'*. As shown 
by Chang^^ and exploited in Sec. IIIII under certain con- 
ditions Lowdin orbitals (LOs) provide the simplest con- 
nection to the k • p method. Another useful property is 
that LOs centered at different sites are orthogonal. 

For the case of a bulk material, SK proceeded to con- 
struct bulk Bloch sums from the LOs 



i.k 



1 



R, 



(A2) 



where N is the number of lattice sites in the crystal and k 
lies within the first Brillouin zone. When the Schrodinger 
equation _ff jV'k) = i'lV'k) is expanded in terms of Bloch 
sums, it becomes an eigenvalue problem 



E 

fi' .i' 



(A3) 



where E^^^ is the energy at point k and we have defined 
a expansion coefficient cf(k) = ('/>fkl^k)- The tight- 
binding matrix elements are easily seen to be given by 



Hi:^ (k) ^ 



E' 



;O|iI|0f,;R, 



(A4) 



In empirical tight-binding (ETB), the matrix elements 
on the right hand side of Eq. (|A4p are taken as adjustable 
parameters 



(R.) 



if;0|i/| 



; Ri 



(A5) 
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The sum over the neighboring sites in Eq. (|A4|I can be 
rewritten as 

(<k 1^1 <k) - E 1^ E ^■'"''E^ (Gr, ) (A6) 

where Yj is the collection of Bravais lattice sites not re- 
lated by point group {Qq) operations G {i.e. sum over 
lattice points belonging to different stars) and Nijj) is 
the number of operations that leave invariant. 

Since G acting on a state |(/>^ \Yj) changes both the 
type of state and its lattice point, it is easy to see that 
Tcrj = GTr-G~^, where Tr- carries a state ; 0) into 
-jYj). Thus, we have 

E'^^^'{Gr,)^(^f;0\HTGr,\(b'::;0) = 

(0f;O|iJGrr^G^i|0^';o). (A7) 

Since G G ^o, we have that HG = GH, therefore we 
arrive at 

E^'iGv,) = D{G)lE^:!::{,v,)D{G~yf,,, (A8) 

where we make use of 

G|0f;r,) = |0^;Gr,)7^(G)^, (A9) 

and D(G)^^ is the element for the matrix corresponding 
to operation G of the fi-th. irrep of the crystal point group. 
We substitute Eq. (|A8|I into Eq. (|A6() to obtain 

rj G 

(AlO) 

where D{Gy^'^ is the representation of G for polar vec- 
tors, which in general will not be an irrep of the point 
group. This form allows for an easier computer imple- 
mentation since there is no need to keep track of which 
neighbor positions have been visited and which not, and 
a single loop over all symmetry operations can be used. 

SK present in their tables results for the Hamiltonian 
matrix elements computed with a method analogous to 
(|A10p . but only single group irreps are used, which ef- 
fectively does away with spin. In Sec. IIIII we carry out 
calculations with double group irreps, thus taking into 
account all spin effects within the single particle approx- 
imation. 

Note that this procedure can deal with general kinds 
of crystal lattices, not only those with a single atom per 
primitive cell, because any symmetry reduction due to 
the atom basis attached to each lattice site is mimicked 
by the imposed symmetry of the Lowdin orbitals. For 
example, in a zinc blende structure orbitals transforming 
according to Td (i.e. without definite parity) are attached 
to fee lattice sites. 



The inclusion of more than one type of orbital per 
primitive cell [i.e. need to have orbitals attached to 
atomic sites rather than lattice sites) is easily achieved 
by introducing and extra index a for the atom type in 
the states and having Rj Rj + d^/ — in Eq. I|A4|I . 
where do, is the offset of species a with respect to the 
lattice site. After that the derivation until Eqs. IjASp and 
(|A10p continues unchanged. 

APPENDIX B: SYMMETRY CONSIDERATIONS 

It is useful to use all the available symmetries to re- 
duce as much as possible the number of independent 

-^ss' This becomes specially important when ETB 

is applied to lower dimensionality structures, as then the 
sum (|A4|) is not carried out for some directions (i.e. pla- 
nar, linear or point orbitals are constructed). 

a. Point group symmetry 

When a operation G is such that Gvj = Vj for a 
given Tj , we obtain a number of consistency requirements 

for the neighbor matrix elements E'^f (r j ) from Eq. HA8|I : 
E^(r,) = DiG)'^^E^^^:{r,)DiG-'rJ,. (Bl) 

b. Time reversal 

For double group irreps, the basis states can be la- 
beled24 by i = —q^, . . . , — 1/2, 1/2, . . .^q^. For example, 
the top of the valence band in zinc blendes would have 
q^ = 3/2, the bottom of the conduction band — 1/2, 
etc. We employ the phase conventions^ 

e|0r;r,)^|,^r;r,) = (-1)^^-10^,; r,) (B2) 

for the action of the time reversal operator 8. Since 
is antiunitary^^, we will have 

E^iji) = (-l)'-+V(-l)*+^'iJ^(_/(r,.)- (B3) 

c. Rotation followed by translation 

If there is an operation such that Grj = — r^, we can 
combine Eq. IjASp with the translational symmetry re- 
quirement E'j^f (— Tj) = -Eiv j^ (rj) to obtain 

El^^*(r,) = i?(G)r,<';'(r,)i?(G-i)';;. (B4) 

Now we can remove the complex conjugation operation 
by employing time reversal, Eq. (jBSp . The result is 

(_l)..+V(_i)»+'i?(G)f,<';'(r,)i?(G-i)^,;. (B5) 
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If the crystal has an inversion center, we will always 
be able to pick G = I, where / is the inversion operator. 
Then, since every irrep /i of a crystal with inversion will 
have a definite parity tt^ , we arrive at 

E'tl_,(r,) = i-iy^+'^^^'i-iy+^'7:,n,,E^^'ir,). (B6) 

d. Block sums 

It is also convenient to state explicitly the transforma- 
tion properties of Bloch sums and their matrix elements. 
This can reduce the number of elements that need to be 
specified. Combining Eqs. IjAip and l|A9|l it is easy to see 
that 

G|<u> = l<Gk>^(G')^., (B7) 
which leads to the relationship^ 

H^(k) ^ D{G)>t,M^: {G-'\.)D{G-')i„ (B8) 

constraining the values of the Hamiltonian matrix ele- 
ments between Bloch sums. 



APPENDIX C: MODELS WITH A FINITE 
NUMBER OF BANDS 

The procedure described in Appendix ^ will only 
yield the exact Hamiltonian when an infinite number of 



Lowdin orbitals (LOs) is included and interaction with 
all neighbors is accounted for. Of course, the localized 
character of the LOs will limit the neighbor distance at 
which there is non-negligible overlap. On the other hand, 
since the number of included neighbors is directly related 
to the number of Fourier components of the bands, the 
effect of coupling to bands outside the model — which will 
add to the number of Fourier components — can be mod- 
eled by increasing the range of the effective Hamiltonian 
interaction^. 

In the extreme case where only two bands are included 
(say, spin up and spin down) , a large number of included 
neighbors will provide a good description of the bands. 
In this context Kramers degeneracy is easily shown for 
doubly degenerate irreps (i.e. doubly degenerate at the 
zone center, but not a priori at a general k point). 
Labeling /i — s- 2q^, time reversal combined with inver- 
sion symmetry [Eq. (|B6(l ] implies £'||(rj) = i?||(rj) and 
£;|i(rj) — E^l{vj) = 0. Thus, the spin 2x2 Hamiltonian 
takes on the form 



n _ /ifTT(k) \ 
^inv-^^ i?TT(k)y' ^ ' 



which shows that the bands will indeed be degenerate at 
a general k point. 



